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We propose to design multispin quantum gates in which the input and output two-state 
systems (spins) are not necessarily identical. We describe the motivations for such studies 
and then derive an explicit general two-spin interaction Hamiltonian which accomplishes 
the quantum XOR gate function for a system of three spins: two input and one output. 
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1. Introduction 

Dimensions of semiconductor computer components will soon reach^ about 0.25 /um 
= 2500 A, still well above the sizes at which quantum-mechanical effects are im- 
portant. However, it is generally expected that as the miniaturization continues, 
atomic dimensions will be reached, perhaps, with technology different from today's 
semiconductors. The physics of quantum-mechanical computation has attracted 
much attention recently.^ Quantum computer is a quantum-coherent system that 
functions as a programmable calculational apparatus, quite unlike its classical coun- 
terparts. It can perform certain tasks^ much faster than the classical computer: the 
quantum interference property yields^ the fast-factoring (Shor's), as well as certain 
other fast algorithms. Recent theoretical results have included identification of 
universal reversible two-bit gates^ and advances in error correction.^ There have 
also been experiments and proposals for experiments,^ realizing the simplest gates. 
However, the actual construction of a macroscopic computer out of a large num- 
ber of quantum bits (two-state systems, qubits) is ellusive^ at the present stage 
of technology. The main obstacle is the sensitivity of coherent quantum evolution 
and interference to undesirable external interactions such as noise or other fail- 
ures in operation, ^'^'^ even though a number of error correction schemes have been 
proposed.^ 

We therefore propose an alternative approach along the lines of the "classical" 
analog computer design, of operating the computer as a single unit performing in 
one shot a complex logical task instead of a network of simple gates each performing 
a simple "universal-set" logical function. The computer as a whole will still be 
subject to errors. However, these will not be magnified by proliferation of sub-steps 
each of which must be exactly controlled. 

Indeed, quantum (and more generally reversible) computation must be exter- 
nally timed: the time scale of the; operation of each gate is determined by the 
interactions rather than by the relaxation processes as in the ordinary computer. 
Furthermore, gate interactions must be externally switched^ on and ofl: because the 
gates affect each other's operation. Time dependence smoother than the on/off 
protocol is possible: see Ref. 8. 

In fact, we consider it likely that technological advances might first allow design 
and manufacturing of limited size units, based on several tens of atomic two-level 
systems, operating in a coherent fashion over sufficiently large time interval to 
function as parts of a larger classical (dissipative) computer. We would like these 
to function as single units rather than being composed of many gates. 

While in principle in a reversible computational unit input and output spins 
(qubits) need not be different, for larger units interacting with the external world 
it may be practically useful to consider input and output separate (or at least not 
identical). Indeed, the interactions that feed in the input need not necessarily be 
identical to those interactions/measurements that read off the output. 

We consider in this work a spatially extended XOR gate based on three spins: 
two input and one output. Generally, we have to address a complicated set of 
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problems: can multispin computational units bo designed with short-range, two- 
particle interactions? Can they accomplish logical functions with interactions of 
the form familiar in condensed-matter or other experimental systems? These and 
similar questions can only be answered by multispin-unit calculations which will 
have to be numerical. Analytical results arc limited to the simplest gates such 
as NOT and XOR, the former studied in Ref. 8, and they provide only a partial 
picture. 

We emphasize that the XOR fmiction is used only as a solvable example of 
a gate with more than two spins, in which one can seek to accomplish a useful 
logical function solely with two-spin interactions. The XOR function can be also 
realized with two spins (one of the inputs serving as the output), for instance as 
the sub-result of the controlled-NOT gatc.^"'' '' It is also important to emphasize 
that while moving from two spins to three spins brings in the issue of the two-body 
interactions, the other important aspect of multispin gates: having the interactions 
short-range, can only be explored with larger systems. 

In Section 2, we define the problem and introduce some notation. In Section 
3, we analyze the matrix forms of the unitary evolution operator and Hamiltonian 
operator. The latter is explicitly calculated in Section 4 and then further refined in 
Sections 4 and 5 to yield a two-spin- interaction result. 

2. The Three-Spin XOR Gate 

We will use the term "spin" to describe a two-state system, and we will represent 
spin- 5 -particle spin-components (measured in units ofh/2) by the standard Pauli 
matrices ax,y,z- We denote by A, B, C the three two-state systems, i.e., three spins, 
involved. We are particularly interested in such initial conditions, at time t, that 
the input spins A and B are in one of the basis states \AB) = |10), |01), or 
1 00), where 1 and refer to the eigenstatos of the z-component of the spin operator, 
with 1 referring to the "up" state and referring to the "down" state. We use this 
convention for consistency with the classical "bit" notion. The initial state of C is 
arbitrary. 

We would like to have a quantum evolution which, provided A and B are initially 
in those basis states, mimics the XOR function: 

A B output 
1 1 

10 1 (1) 

1 1 


where the output is at time t + At. As already mentioned, one way to accomplish 
this is to produce the output in A or B, i.e., work with a two-spin system. The 
Hamiltonian for such a system is not unique. Explicit examples can be found in 
Rcfs. 2, 5, 9. In the case of two spins involved, the interactions can be single- and 
two-spin only. 
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An important question is whether multispin systems can produce useful logical 
operations with only two-spin and, for larger systems, short-range interactions. 
Indeed, two-particle short-range interactions are much better studied and accessible 
to experimental probe than multiparticle interactions. As a solvable example that 

addresses the former issue, here wc require that the XOR result be put in C at time 
t + Ai. The final states of A and B, as well as the phase of C are arbitrary. In 
fact, there are many different unitary transformations, U, that correspond to the 
desired evolution in the eight-state space with the basis \ABC) = |111), |110), |101), 
|100), |011), |010), lOOl), |000), which wc will use in this order. The choice of the 
transformation determines what happens when the initial state is a superposition 
of the reference states, what are the phases in the output, etc. 

Let us consider first, for illustration, the following Hamiltonian, which elludes 
to our more general results below, 



It is written here in terms of the spin components; the subscripts A, B, C denote 
the spins. In the eight-state basis specified earlier, its matrix can be obtained by 
direct product of the Pauli matrices and unit 2x2 matrices I. For instance, the 
first interaction term is proportional to 



etc. This Hamiltonian involves only two-spin-component interactions. In fact, in 
this particular example A and C only interact with B. 

One can show that the Hamiltonian (2) corresponds to the XOR result in C at t+ 
At provided A and B where in one of the superpositions of the appropriate "binary" 
states at t (we refer to superposition here because C is arbitrary at t). There are two 
ways to verify this. Firstly, one can diagonalize H and then calculate the unitary 
evolution operator (matrix) U in the diagonal representation by using the general 
relation (valid for Hamiltonians which arc constant during the time interval At; see 
Ref. 8 for a formulation that introduces a multiplicative time dependence in H), 



and then reverse the diagonalizing transformation. The calculation is quite cum- 
bersome. 

The second, more general approach adopted here is to "design" a whole family 
of two-spin-interaction Hamiltonians of which the form (2) is but a special case, by 
analyzing generally a family of 8 x 8 unitary matrices corresponding to the XOR 
evolution. This "design" program is carried out in the following sections. 

3. The Structure of the Unitary Matrix and Heimiltonian 

We require any linear combination of the states and |110) to evolve into a 

linear combination of 1 1 10) , |100), |010), and |000), compare the underlined quantum 





O-zA ® CTyB <Xi Ic 



(3) 



U = exp {-iHAt/h) 



(4) 
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numbers with the first entry in (1), with similar rules for the other three entries in 
(!)• 

In the matrix notation, and in the standard basis introduced earlier, namely, 
\ABC) = |111), |110), |101), |100), |011), |010), |001), |000), the most general XOR 
evolution operator corresponding to the Boolean function (1), with the output in 
C, is, therefore 
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(5) 



The condition of unitarity, UV^ = 1, reduces the number of independent param- 
eters. Still they are too numerous for the problem to be manageable analytically; 
recall that each nonzero element Ukn is complex and therefore involves two real 
parameters. Thus, we are going to consider a subset of operators of the form (5). 

From our earlier work* we know that one convenient way to reduce the number 
of parameters and at the same time ensure unitarity is to have a single phase factor 
in each column and row of the matrix. Furthermore, we choose a form which is 
diagonal in the states of the A-spin, 



U = 



^4x4 
04x4 



04x4 
W4x4 



(6) 



Thus, A and B arc not treated symmetrically. Here 04x4 denotes the 4x4 matrix 
of zeros. The 4x4 matrices V and W are parametrized as follows: 



V = 



W 
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(7) 



(8) 



The reasons for this choice of an 8-parameter unitary matrix U will become 

apparent in the course of the calculation. Some of the features can be explained 
at this stage as follows. We note that, omitting the direct-product symbols and 
replacing unit matrices by 1, etc., the matrix U in (6) has the structure 



2U={1 + a,A) V + {1- (j,a) W = V + W + <j,a{V - W) , (9) 
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where V and W are operators in the space of B and C . Since U was chosen diagonal 
in the space of A, the Hamiltonian H will have a similar structure, 

2H = P + Q + a,A{P-Q), (10) 

with the appropriate [B C)-space Hamiltonians P and Q. Now in order to avoid 

three-spin interactions, P — Q must be linear in Pauli matrices. On the other hand, 
we also prefer to avoid single-spin (external-field) interactions. Thus, P + Q must 
contain only terms of the second order in the spin components while P — Q must 
contain only terms of the first order in the spin components. This suggests avoiding 
putting phase factors on the diagonal, which would lead to matrices similar to those 
encountered in the NOT-gate calculations^ that are known to be of a structure 
undesirable here: they contain a mixture of first-order and second-order terms. 
The off-diagonal choices remaining are considerably limited; the forms (7) and (8) 
are thus nearly unique. 

In summary, while the arguments are admittedly vague and they do involve 
a certain level of guess, trial and error, the presented parametrization offers a 
good chance that with further restrictions on the parameters a two-spin interaction 
Hamiltonian can be obtained. As will be seen later, five conditions are imposed so 
that the resulting Hamiltonian depends on three (real) parameters. 



4. The Hamiltonian Matrix 

Let us define 



a + P + -f + 5 
M= ^ , (11) 



_ p + io + £,+ri 
and also introduce the reduced operators p and q according to 



(12) 



P = -^^V and g = -Ag. (13) 



Then (4) reduces to 



V = exp(ip) and W = exp{iq) . (14) 

The following calculations are rather cumbersome. Only the results will be 

presented. The algebraic steps omitted are straightforward. First, we diagonalize 
V and W: we calculate their eigenvalues and also the matrices of their normalized 
eigenvectors. The latter can be used to transform to the diagonal representations. 

Specifically, the eigenvalues of V are e''^, ie''^ — e*'', —ie^^. The appropriate 
eigenvalues of p then follow form (14) as /z, -I- ^tt, /i + tt, /i + |7r. Arbitrary mul- 
tiples of 27r can be added to these choices. However, there are certain nonrigorous 
arguments^ for generally keeping the spread of eigenvalues of the Hamiltonian as 
small as possible. Thus, we choose the simplest expressions. The eigenvalues of q 
are determined identically, with ji replaced by v throughout. 
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The next step is to apply the inverse of the diagonahzing transformations for V 
and W to the diagonal 4x4 matrices for, respectively, p and q. The latter contain 
the eigenvalues of p and q as the diagonal elements. The results are the matrix 
forms of the operators p and q in the original representation: 



4 

-p = 



-(l-z)e*("-'^) 
-(1 +i)e*(''-*) 



-(H-i)e*(''-«) 

_e'(/3+7-2/:i) 

-(1 - i)e'(T-'') 



(1 - i)e^^^-^'^ 
Im + 3 



_gi(2AJ-/3-7) 



-(l + i)e' 



i(/:i-/3) 



_(l_i)ei(/3-M) 



fM + 3 / 



(15) 



-9 

TT 



-(l + i)e'('^-'') 
_(l_i)ei(«— ) 

y _gj(2i^-^)-a;) 



_gi(2i/-w-77) 

-(l + «)e*('^-'^) 



-(l + i)e^(''-€) 

_gi(w+T;-2i/) 

^ + 3 
-(1 - j)e'(''-'^) 



_gi(p+w-2i/) 

-(1 - i)e'('^-'') 
-(1 + i)e*('^-'') 



(16) 



5. The Two-Spin-Interaction XOR Hamiltonian 

Thus far we decreased the number of independent parameters in the general unitary 
transformation and chose it to be diagonal in the vl-space. We now "refine" our 
design of the Hamiltonian to favor interactions of the second order in the Pauli 
matrices. First, we note that both P and Q are constant-diagonal matrices. In 
terms of the Pauli matrices, then, both their sum and difference in (10) will involve 
constant terms. These are undesirable because in <7za{P — Q) they lead to terms of 
order one (instead of the desired two), in H, while in P + Q they load to an additive 
constant in H which only affects the overall phase of the unitary transformation 
and is of no interest otherwise. Therefore, we put 



(17) 



in order to nullify these diagonal elements in both P and Q. 

Let us now focus our attention on P — Q which, by (17), is now a matrix with 
zero diagonal. We must impose conditions on the parameters to make P — Q of order 
exactly one in the Pauli matrices. We note, however, that due to zero-diagonal, it 
cannot contain terms. The general form linear in (Jx,cry is 



P-Q=Ie 



X 
X* 
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Y* 
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X* 

Y* 
\ 



X 
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Y 





X 


X* 






, (18) 



where the stars denote complex conjugation, X and Y are arbitrary (complex) 
numbers, and I stands for the unit matrix as before. Thus we require that P — Q be 
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of the form suggested by (18). This imposes several rather cumbersome algebraic 
conditions: two above-diagonal elements of the difference must be equal to zero 
while the remaining four elements must be equal pairwise. After a lengthy but 
straightforward algebra not reproduced here, we conclude that these conditions can 
be satisfied if oi.,0., 7 arc kept as three independent (real) parameters while the 
remaining angles are given by 



5 = — Stt — q; — /? — 7 , 



p = -7r + /3, 



LO = —2-K — a — /? — 7 , 



(19) 
(20) 
(21) 



^ = -7r + 7, 



(22) 



= TT + a . 



(23) 



These conditions take care of the form of P — Q. Interestingly, our results 
below also show that P + Q contains only two-spin interactions with this choice 
of parameters. We have no simple explanation of this property (of the absence of 
first-order terms in P -|- Q). It is probably related to the fact that the structure 
pattern of the original matrices V and W is quite similar even though the precise 
positioning of nonzero elements in them is different. Note that (17) is built into 
(19)-(23). The explicit expressions are obtained by a lengthy calculation, 
P + Q = -^^x 



I 











gi7 _ gi(a-|-/3-|-7) \f2e^^^^^^ 



gi(a-|-/3-|-7) _ gi7 



P-Q = - 



\/27rhi . 



4At 
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gia ^ g-i/3 



(24) 
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(25) 

Finally, wc expand these matrices in terms of products of the Pauli matrices and 
collect terms according to (10) to get 
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H = - ^|\/2 (sin a + sin/3) 0-^AO-a;C - \/2 (cos a - cos/3) cr^ACTj/c 

+V2[smj + sin(a + /? + j)]azAO-xB - a/2 [cos 7 + cos(a + /3 + ^)]azAO-yB 

+\/2 (sin a - sin /3) a^BO-xC - \/2 (cos a + cos /3) azso-yc 

-a/2 [sin 7 - sin(a + /3 + 7)]CTa;Bcr^c + V2[cos7- cos(a + /3 + 7)] o-j^bct^c 

- [sin(a + 7) + sin(/J + 7)] axB<JxC + [cos(a + 7) - cos(/3 + 7)] (JxB<JyC 

+ [ cos{a + 7) + cos(/3 + 7)] ayBO-xC + [ sm{a + 7) - sin(/3 + 7)] ayBCfyC | • (26) 

Note that (2) corresponds to the parameter choice a = /3 = 7 = 0. The Hamil- 
tonian (26) describes the three-spin XOR for arbitrary parameter values. All the 
interactions involved are two-spin as desired. 

The result, however, is not symmetric in any obvious way. It seems to corre- 
spond to complicated tensor interactions involving expressions of order two in the 
components of the three spins involved. No rotational or other symmetry in the 
three-component spin space, or planar symmetry, or uniaxial coupling, arc appar- 
ent. These would correspond to the familiar Hciscnberg, XY, and Ising couplings 
in condensed matter physics. Thus, in order to realize interaction (26) in mate- 
rials, a rather anisotropic medium with highly nonsymmetric tensorial magnetic 
interactions will be required. 

In this respect our analytical attempt to "design" a multispin quantum gate in 
this work may indicate that different roots to the derivation of Hamiltonians should 
be also explored. One could start with the more conventional magnetic interactions, 
isotropic (Hciscnberg) , planar (XY) . uniaxial (Ising) , write down general-parameter 
Hamiltonians, and then adjust the coupling parameters numerically in search of 
those values for which useful Boolean gate operations are carried out. There is no 
guarantee that such a program will succeed. We intend to pursue both approaches 
in our future work. 

In summary, we derived a three-parameter family of Hamiltonians that cor- 
respond to the three-spin XOR gate. While our calculation demonstrates that 
multispin gates can accomplish quantum-logic operations with two-particle inter- 
actions, our results seem also to call for further work seeking improvement in two 
ways. Firstly, our derivation is not general and it has involved a good deal of guess 
work. Secondly, the terms in the resulting Hamiltonians have no obvious grouping 
by symmetries. 
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